Simple analytic model of wormhole formation 



Hidcki Macda* 

Centra de Estudios Cientificos ( C'ECS), Arturo Prat 514, Valdivia, Chile 
(Dated: February 2, 2009) 

An explicit and simple solution representing the wormhole formation is presented. The spacetime 
is constructed by gluing the Minkowski and Roberts spacetimes at null hypersurfaces in a regular 
manner. The parameters in the Roberts solution are required to give the negative kinetic term for 
the massless scalar field. Although a curvature singularity appears at the moment of the wormhole 
formation, it disappears instantaneously. This instantaneous singularity is weak in the senses of 
both Tipler and Krolak along radial causal geodesies. 
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I. INTRODUCTION 

Along with black holes, wormholes are intriguing ob- 
jects in general relativity which have been attracting peo- 
ple even not working in gravitational physics. A worm- 
hole is locally characterized by a "throat," i.e., a two- 
dimensional compact spatial surface of minimal area on 
an achronal hypersurface, connecting some asymptotic 
regions or infinities. Wormholes admit the (apparent) 
superluminal travel as a global effect of the spacetime 
topology [1-3]. Moreover, they are available to make 
time machines [4, 5]. (The readers should refer to [2] for 
a standard textbook and [3] for a nice recent review.) 

The Morris-Thorne static traversable wormhole con- 
necting two asymptotically flat spacetimes is now a well- 
known classic in general relativity [6]. (Static worm- 
hole metrics were obtained even before Morris and 
Thorne [7].) It is known that an exotic matter vio- 
lating the null energy condition is necessary for static 
traversable wormholes in general relativity [2, 8, 9]. This 
is also a natural consequence of the topological censor- 
ship in the asymptotically flat case [10]. Thus, to con- 
struct wormhole solutions with small or even without vio- 
lation of the energy condition has been a big challenge in 
wormhole physics [11, 12]. In fact, it was shown that the 
wormhole spacetime can be constructed with an arbitrar- 
ily small amount of matter which violates the averaged 
null energy condition [13]. This result suggests that the 
wormhole configuration could be realized in the universe 
by some quantum effects violating the energy conditions. 

Then, a natural question is the stability of the worm- 
hole solutions. A wormhole could be formed from the 
gravitational collapse of matter fields possibly violating 
the energy conditions. Also, it may be formed by some 
quantum tunneling effect. The stability analysis is im- 
portant in order to clarify the stable stationary config- 
uration of a wormhole. In the case of the static worm- 
hole solution with thin shells, there exist linearly stable 
configurations depending on the parameter(s) of the so- 
lution [14]. On the other hand, no stable and analytic 
wormhole solutions have been reported in the studies of 
mode analyses and numerical simulations so far [15-18]. 



Independent of the stability of the stationary worm- 
hole configurations, the formation of a wormhole is a 
highly nontrivial problem because it is a dynamical pro- 
cess of the topology change. Actually, the dynamical 
aspects of wormholes have not been well understood so 
far. Although there is a lot of static wormhole solutions 
obtained in the literature, there are few works on the ex- 
act model of wormhole formation from the regular initial 
data. Because the formation or the growth of a wormhole 
is essentially a quite complicated dynamical and inhomo- 
gencous process, numerical methods have been often used 
to study such problems. In these surroundings, exact an- 
alytic models are important to give a transparent picture 
of the phenomenon. They become test beds for the fu- 
ture research and should be intensively investigated to 
complement the numerical works. 

The purpose of the present paper is to give a simple 
analytic model of wormhole formation with a massless 
ghost scalar field. The rest of the present paper is con- 
stituted as follows. In the following section, basic equa- 
tions and a review of the Roberts solution are presented. 
In Sec. Ill, we construct our model and study its global 
structure. In Sec. IV, properties of the curvature singu- 
larity in the Roberts spacetime are studied. Concluding 
remarks and discussions including future prospects are 
summarized in Sec. IV. In Appendix A, the relation be- 
tween the Roberts and the Gutman-Bespal'ko solutions 
is explicitly shown. In Appendix B, the global structure 
of the Roberts solution for the nonghost case is reviewed. 
We adopt the units such that c = G ^ 1 . The metric sig- 
nature convention is taken to be (— , +, +, +), and greek 
indices run over all spacetime indices. The conventions 
of the curvature tensors are [Vp, Wcr]V^ = TZ^iypaV^ and 
V — VP 



II. MODEL AND THE SOLUTION 



We begin with the following action: 



S = I d^x^J—g 



(2.1) 
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where e = 1 and —1 respectively correspond to the real 
and ghost massless scalar field. The energy-momentum 



2 



tensor for a scalar field is given by 

The Einstein equation is 

while the equation of motion for (p is 

ns = 0. 



(2.2) 



(2.3) 



(2.4) 



We see from the basic Eqs. (2.3) and (2.4) that if the 
scalar field in one solution with e = 1 is purely imaginary, 
it can be interpreted as a solution with a ghost scalar field 
(e = -l). 

In this paper, we consider the spherically symmetric 
spacetime (A^^, g^^i,) which is a warped product of a two- 
dimensional constant curvature space (S*^, jij) and a two- 
dimensional orbit spacetime {M'^,gAB) under the isome- 
try of {S'^jjij). Namely, the line element is given by 



ds^ = gABdx'^dx" + Rix'^fdn 



(2.5) 



where A,B^ 0,1; i,j = 2,3 and dfl^ := •jijdx'^dx^ = 
d9^ + s\T? Odip^ . Here i? is a scalar on {M^,gAB) with 
i? = defining its boundary, and 7^ is the unit metric on 
(S'^,7y ) with its sectional curvature fc = 1. The Misner- 
Sharp mass [19] is defined by 



'TIMS 



R 



(l-i?,Ai?'^). 



(2.6) 



Under the assumption that {M'^,gAB) is Minkowski, 
the general homothctic self-similar spherically symmetric 
solution for a masslcss scalar field with e = 1 is given by 

ds^ = -^2dudv + {-uv + CW + C2u'^)dn^, (2.7) 

where Ci and C2 are real constants. For C1C2 = 1/4, it 
is the Minkowski spacetime. For C1C2 7^ 1/4, the scalar 
field is given by 



± 



1 



20F 



arctanh ( 



2Ci{v/u) 



V x/1 - AC1C2 



for Ci^O and 



± 



4V^ 



In 



Co 



(2.8) 



(2.9) 



for Ci = 0, where the value of the constant 0o is mean- 
ingless. The Misner-Sharp mass (2.6) is given by 



rriMS 



(1 - ACiC2)uv 



The Kretschmann invariant K is given by 

(-W + ClU^ +C2m2)4- 



(2.10) 



(2.11) 



The expression (2.8) is convenient to understand the 
codomain of </) because arctanh(zi;) is real, complex, and 
purely imaginary for < (Re w)"^ < 1 with Im w = 0, 
(Re w)^ > 1 with Im w = 0, and Re w ~ 0, respectively. 
The scalar field is real for 



[l-2C,{v/u)r 

- 1 - ACIC2 



(2.12) 



in the case of C1C2 < 1/4 otherwise </> becomes complex. 
We write the condition (2.12) in terms of the square of 
the areal radius = —uv + Civ^ + C2U^ as 



< 



4Cii?2 



u2(l-4CiC2) 



1 < 1. 



(2.13) 



Therefore, under C1C2 < 1/4, the scalar field is real and 
complex for Ci < and Ci > 0, respectively. In the case 
of C1C2 > 1/4, on the other hand, Eq. (2.8) with 00 = 
can be rewritten as 



±- — ;= arctan 



■2V^ 



2Ci (!>/-») 
V4CiC2 - 1 



(2.14) 



where := —1. Therefore, the scalar field is purely imag- 
inary corresponding to a ghost scalar field for C1C2 > 
1/4. 

Here we must give some comments on the history of 
this solution. The solution with Ci 7^ was obtained by 
Roberts in 1989 [20]. Unfortunately, the metric in the 
double null coordinates was erroneously written and the 
correct form was found later by several authors [21-24]. 
In this case, we can set |Ci| = 1 without loss of general- 
ity by the coordinate transformations v := -^Z | Ci | u and 
u := u/ \/\Ci\, so it is a one-parameter family of solu- 
tions. On the other hand, the solution with Ci = was 
obtained by Brady in 1994 [23]. (See also [25, 26].) In 
fact, we can show that the metric of the solution found by 
Gutman and Bespal'ko for a stiff fiuid in 1967 [27] covers 
half of the spacetime (2.7). (See Appendix A.) This is 
because a massless scalar field is equivalent to a stiff fluid 
if the gradient of the scalar field is timelike [28]. Keep- 
ing in mind the history, we call this solution the Roberts 
solution in the present paper. 



III. AN ANALYTIC MODEL OF WORMHOLE 
FORMATION 

In this section, we construct a simple analytic model 
of wormhole formation by gluing the Roberts and 
Minkowski spacetimes in a regular manner. We focus on 
the case of C1C2 > 1/4 corresponding to a ghost scalar 
field, which is required for this construction. The prop- 
erties of the Roberts solution in the case of C1C2 < 1/4 
is reviewed in Appendix B. 

In the case of C1C2 > 1/4, Ci > and C2 > are re- 
quired for the areal radius to be nonnegative. The areal 
radius becomes zero only at u — v = 0. Thus, it is seen 
in Eq. (2.11) that only u ~ v ~ may be a curvature 
singularity. We also sec in Eq. (2.10) that the region with 
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uv > (<)0 has positive (negative) mass. On a null hy- 
persurface of m = or i; = 0, the Kretschmann invariant 
and the quasilocal mass are zero and the derivative of the 
scalar field becomes null. 

The trapped region is given by i? < 2mMS- Since 
"IMS < is satisfied for uv < 0, the trapped region is lo- 
cated in the region of > 0. The trapping horizon [29] 
defined by i? = 2toms is given by 



basis vectors of S defined by := dx^/dy"' arc given by 



4Ci ' 



(3.1) 



for Ci ^ 0, while u = and v = 2C2U for Ci = 0. 
Thus, there are two timelike trapping horizons (3.1) for 
C1C2 > 1/4 with Ci, C2 > in the region of uu > 0. 

It is easy to know the global structure of the Roberts 
solution because {Al'^,gAB) is the Minkowski spacetime. 
u and V are affinc parameters along radial null geodesies, 
and then null infinities are represented by u ±00 or 
V —> ±00. The Penrose diagram of the Roberts spacetime 
for C1C2 > 1/4 with Ci, C2 > is given in Fig. 1. This 
spacetime represents a dynamical wormhole. (Several 
(quasi-)local definitions of a dynamical wormhole have 
been independently proposed so far [30-32].) 




FIG. 1: The Penrose diagram of the Roberts solution (2.7) 
for C1C2 > 1/4 with Ci > and C2 > 0. A star corresponds 
to an instantaneous curvature singularity at m = w = 0. The 
quasilocal mass is positive (negative) in the regions II and 
IV (I and III). The trapping horizons (3.1) are located in the 
regions II and IV. 



d_ 

dv ' 



'dx^ 



(3.3) 
(3.4) 



The basis is completed by N^dx^ = —dv satisfying 
A^^ej^ = — 1 and iV^ef = on S. The only nonvanishing 
component of the transverse curvature Cab '■— N^^-i^e'^e'^ 
of S is 



(3.5) 



The regular attachment on E requires the continuity 
of hab and Cab on both side of E. Since there is no C2 in 
the expressions of hab and Cab, two Roberts spacetimes 
with the same nonzero Ci but different C2 can be at- 
tached in a regular manner at u = 0. Thus, as a special 
case, the Roberts spacetime (2.7) with Ci = Ci{^ 0) and 
C2 = C2 can be attached to the past Minkowski space- 
time at u = in a regular manner, of which metric is 
given by Eq. (2.7) with Ci = Ci and C2 = l/(4Ci) ^ C2. 
Similarly, it is shown that two Roberts spacetimes with 
the same nonzero C2 but different Ci can be attached in 
a regular manner at 1; = 0. 

By gluing the Roberts spacetime with 4C1C2 > 1, 
Ci,C2 > and the Minkowski spacetime(s) at u — 
and/or w = in a regular manner, we can construct 
spacetimes representing wormhole formation from the 
initial data with a regular center. The Penrose diagrams 
for these spacetimes are given in Fig. 2. 

The attachment of the Roberts spacetime to the 
Minkowski spacetime in the case of the ghost scalar field 
has been mentioned in [35] without detailed calculations. 
It is claimed there that the instantaneous singularity 
u = V = in the Roberts spacetime can be removed by 
gluing the Minkowski spacetime at m = or w = 0. Ob- 
viously, the curvature invariants do not blow up if an ob- 
server approaches there from the Minkowski region, how- 
ever, they certainly blow up along some causal geodesies 
emanating from u = w = in the Roberts region. As a 
result, there is still a naked singularity at m = w = in 
the resulting spacetime. The details will be presented in 
the next section. 



IV. PROPERTIES OF THE INSTANTANEOUS 
SINGULARITY 



Now we show that the Roberts spacetime can be at- 
tached to the Minkowski spacetime at m = or u = in a 
regular manner, i.e., without a massive thin shell on the 
hypersurface. (See [33, 34] for the matching condition 
on a null hypersurface.) We consider a null hypersurface 
M = as a matching surface, which we call E. (The ar- 
gument is similar for v ^ 0.) The induced metric hab on 
E is given by 

dsl = habdy^dy^ := Civ'^d^^, (3.2) 

where ?/° = {v,9,ip) is a set of coordinates on E. The 



In the last section, we constructed a spacetime repre- 
senting the wormhole formation. One problem in this 
spacetime is an instantaneous curvature singularity at 
u = V = which appears at the moment of the wormhole 
formation. In this section, we show that it is a naked but 
weak singularity. 

A. Nakedness 

First we show that both radial and nonradial causal 
geodesies emanate from u = v = 0, i.e., it is certainly a 
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FIG. 2; The Penrose diagrams representing wormhole forma- 
tion from the initial data with a regular center. The Roberts 
spacetime with C1C2 > 1/4, Ci > 0, and C2 > (the shad- 
owed region) is attached to the past Minkowski spacetimes at 
(A) u = with i; > and v = with u > 0, (B) u = 0, (C) 
V — with M > 0, and (D) u = with v < and v = with 
u < 0. A thick line corresponds to a symmetric center in a 
Minkowski spacetime. ti represents a spacelike hypersurface 
with a regular symmetric center, while t2 represents a space- 
like hypersurface with distinct spacelike infinities without a 
regular center. 



naked singularity. The Lagrangian to give the geodesic 
equations is 



= -uv + -{-uv + Civ^ + C2u'^){e^ + sin^ Oip^). 

^ (4.1) 

where a dot denotes the derivative with respect to the 
affine parameter A along a geodesic. Because of spherical 
symmetry, we can set = tt/2 without loss of generality. 
The metric (2.7) is independent of (f, so that from the 
Lagrange equation 







d dL dL 



we obtain a conserved quantity along a geodesic as 



(4.2) 



dip 



Then, the geodesic equations (4.2) arc written as 
i-v + 2C2u)$^ 



Q = v + 

= u + 



2(-UU + Ciu2 + (72^2)2^ 

(-u + 2Cit;)$2 

2{-UV + Ciw2 + C2m2)2 



(4.4) 
(4.5) 



The tangent vector of a nonspacelike geodesic := 
satisfies 



(4.6) 



where e is and —1 for null and timclike geodesies, re- 
spectively. This equation is written as 



-2uv + {-uv + Civ^ + C2W^)~^$^ 



(4.7) 



The Roberts spacetime admits a homothctic Killing 
vector (,^{d/dx^) — u{d/du) + v{d/dv) satisfying 



(4.8) 



Then, we obtain 
d 



dX 



S.fj.-.i'k^k 
9fj,uk^k 



(4.9) 



where we used the fact that k^ is tangent to a geodesic 
and Eq. (4.8). Hence we obtain ^^/c^ = Dq + e\, or 
equivalently 



—uv — vu — Dq + eX, 



(4.10) 



{-uv + Civ"^ + C2U^)(p- 



(4.3) 



where Dq is a constant. Equation (4.10) is integrated to 
give 

-™ = Di +D0A+ ieA^ (4.11) 

where Z^i is a constant. 

We are now interested in the geodesies emanating from 
u = V = 0. Without loss of generality, we can set A such 
that A = corresponds to u = tj = 0. Thus, we consider 
the case with Di = 0. Now the geodesic equations reduce 
to 

e = -2uv + {-uv + Ciw^ + C2^t^)~^$^ (4.12) 
-uv = DqX+]^£X^. (4.13) 

We obtain the master equation for u{X) from above equa- 
tions as 

du _ Au{eX + Do) ± u.yA[D'^A + 4A<I>2u2(£A + 2Do)] 
dX ~ AA(eA + 2Dq) ' 

(4.14) 

A ~ 4C2U'' + 2Au2(£A + 2Do) + CiX'^{eX + 2Dof . 

(4.15) 
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First we consider radial geodesies ($ = 0). For the 
null geodesies (e = 0), the solutions of Eqs. (4.12) and 
(4.13) passing through u = v ^ are u = or v = 
0. Along these radial null geodesies, the Kretschmann 
invariant and the quasilocal mass are identically zero. 
For the timelike geodesies (e = —1), the general solution 
of Eq. (4.14) is given by u = uqX and u = uo(2Do — A), 
where uq is a nonzero constant. These two coincide for 
Do = and the latter does not pass through u = for 
Dq 7^ 0. Finally, the solution passing through u ~ v ~ 
is given by 



UqX, 



-A, 



(4.16) 



which corresponds to Dq = in Eq. (4.13). The 
Kretschmann invariant along these radial timelike 



geodesies is given by 

192^8(1 - 4CiC2)^ 



K 



(-2w2 + Ci+4C2u4)4A4' 



(4.17) 



which diverges at A = 0, i.e., u ^ v = 0. Along these 
radial timelike geodesies, the quasilocal mass is given by 



mus 



(1 - ACiC2)uqX 

' 4^-2^2 + Ci+4C2M^^ 



(4.18) 



which converges to zero at A = 0. 

For the nonradial geodesies (<& 7^ 0), there is a solution 
of Eqs. (4.12) and (4.13) passing through u = v = 0, 
which is given by 



-A(2go + eX)[{Dl + 2<1>2) ± ^{D^ + 2<1>2)2 - 4D^^C,C2] 
4D^C2 

-X{2Dq + eX)[{Dl + 2j>2) t ViP^ + 2<i>^)^ - ^DtC,C2] 
^DlCi 



Under 4C1C2 — 1 > and Ci, C2 > 0, the conditions for 
V? and v'^ to be real and positive are Dq <0 and 



4 < ^'^^ - — m — ■ 



(4.21) 



Because Dq and $ are independent parameters which 
characterize a geodesic, the right-hand side of Eq. (4.21) 
varies from 1/4 to infinity. Hence, for any values of Ci 
and C2 satisfying C1C2 > 1/4 and Ci, C2 > 0, there are 
nonradial causal geodesies passing through m = u = 0. 

Along these nonradial geodesies, the Kretschmann in- 
variant diverges around w = w = as 



hm K ~ 37/o^;o(l - AC^C2f 
A^o ~ {-uqvq + Civl + C2ulYX'^ ' 



(4.22) 



-{Dl + 2$2) ^ ^[Dl + 2<^^Y~ADlCiC2 



2DqC2 



(4.23) 



-(Jg + 2$2) ± ^^Dl + 2<1>2)2 - 4/^4^1 C2 
2DqCi 



(4.24) 

while the quasilocal mass converges to zero at u = w = 



(1 - AC,C2)uqVqX^/^ 

hm TTiMS , (4.25) 

^^J-UQVQ + Civl + C2UI 

We note that Eqs. (4.22) and (4.25) are exact expressions 
for nonradial null geodesies. 



B. Strength 



Next we consider the strength of the singularity at u = 
V ^Q. As definitions of the strength of singularities, the 
strong curvature condition (SCC) [36] and the limiting 
focusing condition (LEG) [37] were proposed by Tipler 
and Krolak, respectively. We consider a geodesic (N), 
affinely parametrized by A, with the tangent vector fc^, 
terminating at or emanating from a singularity, where 
A = 0. SCC and LEC imply that N is emanating from or 
terminating in the Tipler's strong and the Krolak's strong 
curvature singularities, respectively [38]. The physical 
content of the Tipler strong is that the volume element of 
physical objects (constructed by the Jacobi fields along 
N) converges to zero at the singularity. On the other 
hand, the physical content of the Krolak strong is that 
the expansion along N diverges at the singularity, but 
still the volume element remains finite. (See also [39, 40] 
for the textbook.) 

The necessary conditions for SCC and LFC are avail- 
able [38]. Let E'fjJl = 1,2,3,4) a parallelly propagat- 



ing frame along N satisfying E^^^E(^i-j^ ~ rj-^^-^rj(^2)ii 



-£^("4)^(4)^ 



= 1 if N is timelike and 



if N is null. All other products vanish and i?(4)p 
If SCC is satisfied along N. then limA^o '^^'^'^''(4)(/)(4) 
does not converge for /, J = 1,2,3 and I, J = 1,2 
in the cases where N is timelike and null, respectively, 
where '^^ ^^'^ Riemann tensor in the paral- 

lelly propagating frame. If LFC is satisfied along N. then 
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lim;^^o AT?. 



(4)(./)(4) 



docs not converge for /, J = 1, 2, 3 



and I, J = 1,2 in the cases where N is timehke and null, 
respectively. 

We show that the singularity at w = u = in the 
Roberts spacetime is weak in the senses of both Tipler 
and Krolak for radial causal geodesies. For the radial 
causal geodesies, of which tangent vector has the form of 
kt" = (fc". A:", 0, 0) =: E^^^ where fc^fc" = is satisfied for 

null geodesies, the angler bases E'^'j-^ {I — 1,2) are given 
as 



E^ — 
E^ — 

(2) dxt- 



RdO' 
1 d 



which satisfy E'^j-^E^^j)^ 



(4.26) 

Sij. The only nonzero compo- 



nent of TZ'' 



(4)1.(4) 



Tl^^^kPk- is 



v\2 



TV 



(4)i(4) 



(4C1C2 - l)(i'fc" - uk") 
uv + Civ^ + C^u^f 
(I) 



A{-uv + 



(4.28) 



Thus, the only nonzero component of TZ (4)(y)(4') 

{AC1C2 - l)(i'fc" - uk^' f 



■r: 



(1) 

(4)(1)(4) 



7^ 



(2) 

(4)(2)(4) 



A{-UV + CiV^ + C2U2)2 ■ 

(4.29) 

These quantities are identically zero both for radial null 
geodesies (m = or w = 0) and radial timelike geodesies 
(4.16). Hence, it is concluded that the singularity at 
M = w = is weak in the senses of both Tipler and 
Krolak for radial causal geodesies. 

For nonradial geodesies, it seems to be cumbersome 



to examine TZ 



(!) 



(4)(J)(4)- 



Instead, we here consider the 



behavior of 7?,(4)(4) TZfj^^k^^k"^ , which is used in the 
sufficient conditions for SCC and LFC. SCC is satis- 



fied if limA^o ^ 



(4) (4) 



> and LFC is satisfied if 



lim^^o A7?.(4)(4) > [38, 40]. The nonzero components 
of the Ricci tensor of the Roberts spacetime are given by 

(4.30) 



uv{4CiC2 ~ 1) 

2{~UV + Civ2 + C2U2)2 ' 

u^{4CiC2 - 1) 
~ 2{-uv + Civ^ + CW)^ 



(4.31) 
(4.32) 



Finally, for nonradial causal geodesies fc^ = (ii,w,0, 1^), 
where u, v, and (p are obtained from (4.19), (4.20), and 
(4.3), respectively, we obtain 7?.(4)(4) = 0, which immedi- 
ately implies A27^(4)(4) = and A7?,(4)(4) = along the 
geodesies. Although this result does not directly mean 
that neither SCC nor LFC is satisfied, it would suggest 
that the singularity at u = = is weak also along 
nonradial causal geodesies. 



V. SUMMARY AND DISCUSSIONS 

In this paper, we constructed an explicit and simple 
model of wormhole formation from the initial data with 



a regular center. The spacetime represents the wormhole 
formation with a massless ghost scalar field. In this con- 
struction, the matter region represented by the Roberts 
spacetime is attached to the past Minkowski spacetimes 
at null hypersurfaces in a regular manner. 

This construction has been mentioned in [35] without 
a detailed analysis. Actually, unlike the authors' claim 
in [35] , a naked curvature singularity appears at the mo- 
ment of the wormhole formation. However, we showed 
that it is only instantaneous and weak in the senses of 
both Tipler and Krolak. This class of weak singularities 
could be harmless because it would be dealt with in some 
distributional sense. 

In this context. Hay ward and Koyama constructed an 
analytic model representing the wormhole "formation" 
from the Schwarzschild black hole [41]. Although their 
model does not contain a singularity at the moment of 
"formation" , it does not represent the wormhole forma- 
tion from the initial data with a regular center, i.e., there 
is no topology change in their model. While they defined 
a wormhole throat by a class of trapping horizons [30, 31], 
there is a wormhole throat on any spacclike hypersurface 
in their model. (See [32] for the discussions of the (quasi- 
)local definition of a wormhole throat on a spacelike hy- 
persurface.) 

In the spatially compact spacetime, the wormhole for- 
mation necessarily requires the appearance of singular- 
ities or closed timelike curves. This is a consequence 
of the result by Geroch about the topology change of 
spacetimes [42]. Hence, under the suitable chronology 
condition, the singularity formation is inevitable. In this 
paper, on the other hand, the wormhole formation in the 
spatially noncompact spacetime is considered. The sin- 
gularity formation would be also inevitable for the worm- 
hole formation in this case, however, as far as the author 
knows, this is an open question. 

The solution presented in this paper will be a simple 
analytic model to study the formation of a wormhole. 
In this context, the stability of the wormhole formation 
is an important future work. While the stability of the 
Roberts solution was studied in the case of the positive 
kinetic term of the scalar field [43], it is still an open 
question for the ghost scalar field. 
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APPENDIX A: THE RELATION BETWEEN THE 
GUTMAN-BESPAL'KO AND THE ROBERTS 
SOLUTIONS 

In 1967, Gutman and Bespal'ko obtained a spherically 
symmetric solution for a stiff fluid, i.e., a perfect fluid 
with an equation of state p = /x, where p and /x are 
the pressure and energy density, respectively [27]. (See 
also [44, 45] for the generalized solution.) The energy- 
momentum tensor for a perfect fluid is given by 



^ P9tii^ + ip + f-i)u^u^, 



(Al) 



where is the four- velocity of the fluid element. The 
Gutman-Bcspal'ko solution is given in the comoving co- 
ordinates as 



-z'^il + Cie^' + C2e~^')dn'^, (A2) 



P 



1 - 4C1C2 



87rz2(i + Cie2*4-C2e-2*)2^ 



(A3) 



where the constants Ci and C2 satisfy C1C2 < 1/4 for 
nonnegative energy density. In the case of 4C1C2 = 1, 
the solution gives the Minkowski solution. The Gutman- 
Bespal'ko spacetime admits a homothetic Killing vector 
^''{d/dx'^) = z{d/dz) satisfying C^g^^ = 2g^^. 

We show that the Gutman-Bespal'ko spacetime covers 
half of the Roberts spacetime. In 1988, Madsen showed 
the equivalence between a massless scalar field and a stiff 
fluid [28]. This is easily generalized to the ghost case as 
shown below. The energy-momentum tensor for a stiff 
fluid is 



(A4) 



If is vorticity free, which is satisfied in the spherically 
symmetric spacetime, one can show that this matter field 
is equivalent to a massless scalar field </>, of which gradient 
(j)^^ is timelikc. The corresponding energy density and 4- 
velocity are given by 



1 



±- 



(A5) 
(A6) 



with which Eq. (A4) coincides with Eq. (2.2), where the 
sign in (A6) is chosen so that is future-directed. 

For the Gutman-Bespal'ko solution, the corresponding 
scalar field with e = 1 is 



± 



1 



2t 



for Ci 7^ and 



2V^ 



± 



arctanh 



1 + 2Cie' 
"VVl-46\6'2 



40F 



In 



(A7) 



(A8) 



for Ci = 0. 

In the Gutman-Bespal'ko spacetime, the metric on 
(A/2^ gj^g) is Minkowski in the Rindler coordinates, while 
it is also Minkowski but in the double null coordinates in 
the Roberts solution. By the coordinate transformations 



t ■ 



T 

arctanh ( — 
A 



(A9) 



of which inverse transformations are 

T = z sinht, X ^ z coslit, 
the two-dimensional Rindler metric ds'^ 



is transformed into dsj - 
null coordinates such as 

T-X 



^dT^ 



2 ' 
dX\ 



V2 ' 



T + X 



(AlO) 

-z'^dt^ + dz^ 
Adopting the 



(All) 



we obtain dsj 



-2dudv. 

Indeed, by the direct transformations 



zk—=ze 
V2 



1 t 



(A12) 



of which inverse is 



t = arctanh 



V + u 



z = ±\/-2uv, (A13) 



the Roberts metric (2.7) is transformed to the Gutman- 
Bespal'ko metric (A2). Therefore, we may call the 
Gutman-Bespal'ko metric (A2) the Rindler chart of the 
Roberts metric. Because of uv < 0, the Rindler chart 
covers only half of the Roberts spacetime. (The regions 
I and III in Figs. 1 and 3.) 

On the other hand, by the coordinate transformations 



T = t coshz, 
of which inverse is 



X ~ t sinh 



i=± Vt2 - A2, S = arctanh ( ^ 



(A14) 



(A15) 



the two-dimensional Minkowski spacetime dsj = —dT^ + 
dX^ is transformed to the Milne form dsj = —dP+Pdz^. 
Thus, by the direct transformations 



u 



1 ~ 
V2 



1 ~ ? 
V2 ' 



of which inverse is 

i = ±V2uv, 



= arctanh 



+ u 

the Roberts solution is transformed to 



(A16) 



(A17) 



ds^ = -dP + Pdz^ 

+ ip(-l + Cie2^~ 



± 



1 



2^/^ 



arctanh 



+ C2e-^'')dn\ (A18) 
l-2Cie2^^ ^ 



(A19) 
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for Ci 7^ 0. For Ci = 0, the scalar field is transformed to 



1 



In 



(A20) 



We may call this metric the Milne chart of the Roberts 
metric. Because of uv > 0, the Milne chart covers the 
regions II and IV in Figs. 1 and 3. This spacetime ad- 
mits a homothetic Killing vector £^^{d/dx^) = t{d/dt) 
satisfying C^g^^ = 2g^^. 



If C1C2 > 0, both of them are timelike or spacelikc, while 
if C1C2 < 0, one is spacelike and the other is timelike. 
For Ci ~ and C2 ^ 0, there are null and nonnull central 
curvature singularities at m = and u = (l/C2)v, respec- 
tively. For C2 = and Ci 7^ 0, there are null and nonnull 
central curvature singularities ai v = and v = (1/Ci)u, 
respectively. For Ci = and C2 = 0, there are null cen- 
tral curvature singularities at u — and v ~ 0. The Pen- 
rose diagram of the Roberts spacetime for C1C2 < 1/4 is 
given in Fig. 3. 



APPENDIX B: THE ROBERTS SOLUTION FOR 

C1C2 < 1/4 

In this appendix, we review the properties of the 
Roberts spacetime with C1C2 < 1/4 corresponding to 
the positive kinetic term of the scalar field. First we see 
in Eq. (2.10) that the region with uv > has negative 
mass. In the case with Ci 7^ 0, C2 7^ 0, and C1C2 7^ 1/4, 
there arc nonnull central curvature singularities located 
at 

- = (Bl) 



As shown in the main text, the Roberts spacetime can 
be attached to the Minkowski spacetime at a null hyper- 
surface u = 0oru = 0ina regular manner if that hy- 
persurface is regular. The resulting spacetime can be a 
model of the gravitational collapse leading to the naked 
singularity formation. This model has been considered 
in the context of critical phenomena or cosmic censor- 
ship [20, 23-25, 43, 46, 47]. 
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